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■ Abstract 

, In an analysis of the IRAS 1.2 Jy redshift catalogue, with emphasis on the separate 

' examination of northern and southern parts (in galactic coordinates), we found that the 

^ clustering of galaxies differs significantly between north and south, showing fl uctuations in the 



^\ . clustering properties at least on scales of 100/i~^Mpc (Kerscher et al. 1997). In Section ^ we 

give a brief description of our morphological method which is based on Minkowski functionals; 
in Section |2| we show the results obtained from the IRAS 1.2 Jy galaxy catalogue. Section |^ 
contains a discussion of several error estimates and in Section ^ we select different subsamples 



O ' from the 1.2 Jy catalogue according to flux and colour and validate the results of Section ^ 

, Furthermore we look closer at the spatial origin of the fluctuation in Section ^ and we compare 

^ ' with optically selected galaxies from the CfAl survey in Section pi 



1 Minkowski functionals 



■ We consider a set {xi} ^^ oi N points give n by the redshift-space coordinates of the galaxies in the 

IRAS 1.2 Jy catalogue ( |Fisher et al. 1995| ). To characterize t he properties of this point set with the 



help of Minkowski functionals introduced into cosmology by Mecke et al. (1994) , we decorate each 
point X j with a ball S, (xi) of radius r and consider the union set AN{r) = Ui=o ^r{^i)- Hadwiger| 



(1957)1 proved that in three-dimensional space the four Minkowski functionals M^=o,i,2,3(^Af(?')) 



give a complete morphological characterization of the body ^Ar(r). For the interpretation of these 



functionals in terms of geometrical and topological quantities see Mecke et al. (1994) and Table 
|l|. Reduced, dimensionless Minkowski functionals $^(^Ar(r)) may be constructed by normalizing 
with the Minkowski functionals Mfj,{Br) of balls, 



For a Poisson process the functionals can be calculated analytically (see Mecke & Wagner 1991) 
with the results: 

$^ = (1 -e-'') rr\ $f = e~'', 

(2) 

$P = c-Ml-^^), <ff = e-" (1-377+ ^^'), 



*Proc. 2"=^ SFB workshop on Astro -particle physics Ringberg 1996, Report SFB 375/P002 (1997), R. Bender, 
T. Buchert, P. Schneider (eds.), in press. 
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Table 1: Minkowski functionals in three-dimensional space expressed in terms of more familiar 
geometric quantities. 

where rj := nMo{Br) = n Aur^ /i, with the mean number density n. As seen from Eqs. (|), the 
Minkowski functionals p,{An {r)) , fJ. — 1,2,3, which are supported by the surface, are proportional 
to e^'''-'"-' for a Poisson process. A similar exponential decay is also found in the case of more general 
cluster processes. Since we are interested in the behavior of the point process on large-scales we 
remove the exponential term and consider functionals 0/^=1.2,3, defined by 

M^N{r)) = ^p^^) ■ (3) 

2 Results from the 100/i^Mpc sample 



We consider a volume limited sample of the IRAS 1.2 Jy galaxy catalogue (Fisher et al. 1995) 
with limiting depth of lOO/i-^Mpc (Hq^IOO h km s^^ Mpc'^). 

Fig. 1^ displays the values of the Minkowski functionals of the southern and northern parts in 
comparison with the functionals of a Poisson process with the same number density. The northern 
part contains 352 galaxies, and the southern part 358 galaxies. The errors of the Poisson process 
were calculated from twenty different realizations. To estimate the error from the catalogue data 
we calculated the Minkowski functionals of twenty subsamples containing 90% of the galaxies, 
randomly chosen from the volume limited subsample (see Section |^). 

In both parts of the 1.2 Jy catalogue the galaxy clustering on scales up to 10/i~^Mpc is dis- 
tinctly stronger than in the case of a Poisson process, as inferred from the lower values of the surface 
functional 01, the integral mean curvature functional (j)2 and the Euler characteristic functional 
(/)3. Moreover, the northern and southern parts differ significantly in their morphological features, 
with the northern part being less clumpy. The most conspicuous features are the enhanced surface 
area 0i in the southern part on scales from 12 to 20 /i~^Mpc and the decrease of the integral mean 
curvature 4>2 which sets in at 14 h~^Mpc. This behaviour indicates that dense substructures in 
the southern part fill up at this scale (i.e. the balls in these substructures overlap strongly); this 
is probably the signature of the Perseus-Pisces supercluster (compare Section ||) . On scales from 
15 to 20 ft.~"'^Mpc, the integral mean curvature (j>2 is negative indicating concave structures. In the 
southern part the Euler characteristic is negative in this range, therefore, the structure is dom- 
inated by interconnected tunnels (negative contributions to the Euler characteristic; completely 
enclosed voids would yield positive contributions). 



3 Error estimates 

To find out how the error from subsampling is related to the intrinsic variance of an ensemble we 
took fixed realizations of a Poisson process within the sample geometry (with the same number of 
galaxies as in the volume limited samples) and calculated the error via subsampling, using again 
90% of the points (see Figure |2|). This error turns out to be about two times smaller than the 
ensemble error estimated from 20 realizations of a Poisson process with the same number density. 
A comparison with the variance of a Poisson process may be helpful, but the fluctuations of the 
functionals calculated for the galaxy distribution exceed this variance (compare the fluctuations 
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Figure 1: Minkowski functionals (/)^ of a volume limited sample with lOOft-^^Mpc depth; dark 
shaded areas: southern part; medium shaded areas: northern part; dotted areas: Poisson process 
with the same number density. The areas correspond to Itr errors, as explained in the text. 
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between north and south with the variance of the Poisson process in Fig. |). With the errors 
calculated from subsampling we only check whether our estimators are robust; the fluctuations 
may be signiflcantly larger. 

To obtain another estimate for measurement errors, we randomized the redshifts using the 
quoted redshift errors as the standard deviation (and using the mean redshift error, if none is 
quoted). These errors are approximately two times smaller than the errors from subsampling. 
Even if we increase the quoted redshift error by a factor of five, the errors in the values of the 
functionals are only of the same order as determined from subsampling (compare Fig. |l| and Fig. |^). 

We also took care for incomplete sky coverage. Apart from a 10 degree wide zone of a voidance. 



whose boundary effects are completely removed according to Bchmalzing et al. (1996) , there are 
additional empty regions in the redshift catalogue due to a lack of sky coverage or confusion 
in the point source catalogue. In the northern part these regions account for 3.2 % of the sky 
in the southern part they account for 4.5 %. To estimate the influence of these regions on the 
statistics we added Poisson distributed points with the same number density. The additional error 
introduced from these random points is much smaller than the errors from subsampling. Moreover, 
no systematic effect is seen, the curves overlap completely. 



4 Selection effects 

Several selection effects might enter into the construction of the catalogue. Therefore, we draw 
subsamples selected according to special features like "colours" or higher limiting flux. 

Before listing these tests, we address the important issue of 'sparse sampling'. This issue is 
especially relevant for the interpretation of our tests, since all these selected samples incorporate 
less galaxies than the volume limited sample with 100/i^^Mpc depth. In Figure ^ we show how 
sparse sampling affects the surface functionals <j)i{AN{'r)) of the volume limited samples with 
100/i~"'^Mpc by only taking a fraction of the galaxies into account (the other functionals behave 
similarly). By reducing the number of galaxies the error increases and the mean values tend 
towards the values for a Poisson process. 

We calculated the Minkowski functionals of a volume limited sample with 100/i~^Mpc depth 
but now with limiting flux equal to 2.0 Jy, with 129 galaxies in the north and 141 galaxies in 
the south. Although the noise increases for larger radii, since fewer galaxies enter, the above 
mentioned features and the differences between northern and southern parts are still detectable 
(Figure I). 

We selected "hot" galaxies, with a flux ratio /loo/Zeo < 1-5, "warm" ones with 1.5 > /loo/Zeo < 
3 and "cold" galaxies with /loo/Zeo > 3 from the 1.2 Jy catalogue; /loo and /go denote the flux 
at 100/i and 60/i respectively. We calculated the Minkowski functionals of volume limited samples 
with lOOft-^^Mpc taken from the "hot" (106 in the north and 116 in the south) and "warm" (239 
in the north and 227 in the south) galaxies (only 7 (north) and 15 (south) galaxies are "cold"). 
Again the error increases, but the features and the differences remain visible (compare Figure |6| 
and 0). With more reflned but essentially similar criteria for distinguishing "warm" and "cool" 



galaxies |Mann et al. (1996) find only a small dependence of clustering on the temperature in the 
QDOT. 



4 





-5 - 



J] 



Figure 2: Minkowski functionals 0^ of Poisson process realizations; the dark shaded area is the 
1(7 error determined with 90% subsaniphng out of one reahzation, the hght shaded area represents 
the 1(7 ensemble error estimated from 20 different realizations. 
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Figure 4: The surface functional (j)i{AN{r)) of a volume limited sample with 100/i~^Mpc depth for 
randomly drawn subsamples with 70% (left figure), 50% (middle figure), and 30% (right figure) 
galaxies. Dark shaded areas: southern part; medium shaded: northern part; dotted: Poisson 
process with the same number density. 
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Figure 5: Minkowski functionals 0^ of a volume limited sample with 100/i^^Mpc depth; now 
the limiting flux is 2.0 Jy. The shaded areas are the la errors determined from randomizing the 
redshifts. 
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Figure 6: Minkowski functionals (pf^ of a volume limited sample with 100/i^^Mpc depth with onb 
"hot" galaxies included (limiting flux 1.2 Jy). The shaded areas are the la errors as in Figure |l| 
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Figure 8: Sketch of the sample cut along the plane through our position, the north galactic pole 
(NGP) and the galactic center (GC). The figure on the right shows how we label them. 



name 


lattitude range 


longitude range 


rn 


5° < < 90° 


0° < I'' < 180° 


In 


5° < fo" < 90° 


180° < < 360° 


rs 


-5° > fo" > -90° 


0° < /" < 180° 


Is 


-5° < &" > -90° 


180° < /" < 360° 



Table 2: The angular ranges covered by the four samples shown in Fig. g|. 

5 Where do fluctuations originate? 

Now we try to find out where the fluctuations have their spatial origin. To do this we cut the north 
an south into two parts each, as shown in Fig. ^ and in Table |2| The four samples approximately 
contain 160 galaxies. 

Fig. ^ shows that all four parts differ from each other. Both southern parts cluster more 
strongly than the northern parts, consistent with Fig. |l|. The strongest clustering is seen in the 
sample Is, the weakest clustering is in the sample In. The Perseus-Pegasus-Pisces Supercluster 
lies in the sample rs. 



6 Comparison with the CfAl 

An objection against results based on IRAS selected samples is that we mainly look at infrared 
active, spiral galaxies. Therefore we compared samples from the CfAl galaxy catalogue with the 
samples from the 1.2 Jy catalogue. We consider volume-limited samples extracted from CfAl 



( pSuchra et al. 1983| ), whose members have apparent magnitudes less than 14.5. The limits of the 
samples considered are galactic latitude 6 > 40° and declination 6 > 0°, the volume limitation is 
performed with 100/i~^Mpc depth leaving 99 galaxies within the sample geometry. The volume 
limited sample of the 1.2 Jy with 100/i~^Mpc has 115 galaxies within the CfAl window. To 
compare both samples we draw random subsamples from the 1.2 Jy to reproduce the number 
density in the CfAl sample. 

The 1.2 Jy galaxy catalogue mainly consists of spiral galaxies, and is therefore undersampling 
the cluster cores. In complete agreement with this, we see in Figure |lO|less clustering in the 1.2 Jy 
on scales up to 10ft.~^Mpc, deduced from the lowered surface functional (pi. On large scales both 
samples are compatible, no systematic difference is seen between the IRAS and the optical sample. 
On small scales both samples differ from the Poisson results but are compatible with Poisson on 
large scales. 

Since we are severely affected by sparse sampling (see Section m we checked the above detected 
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Figure 9: Minkowski functionals (j)^ of a volume limited sample with depth; cut into 

pieces as shown in Fig. ^j. The shading of the la areas (determined from subsampling 90% of the 
galaxies) is getting darker in the following ordering of the samples In, rn, rs, Is. 
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Figure 10: Minkowski functionals 0^ of a volume limited sample with 100/i~"'^Mpc depth of CfAl 
(dark shaded) compared to the Minkowski functionals of 1.2 Jy (medium shaded) and a Poisson 
process (light shaded). The shaded areas are the la errors determined from random subsampling 
90% of the galaxies of the CfAl and 77% of the 1.2 Jy. 
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Figure 11: Minkowski functionals 4>^ of a volume limited sample with 60/i~^Mpc depth of CfAl 
(dark shaded) compared to the Minkowski functionals of 1.2 Jy (medium shaded) and a Poisson 
process (light shaded). The shaded areas are the la errors determined from random subsampling 
42% of the galaxies of the CfAl and 90% of the 1.2 Jy. 



14 



concordance of the morphological properties of optical and IRAS selected galaxies on large scales 
in looking at samples with 60/i~^Mpc depth. Now the CfAl sample includes 215 galaxies and the 
1.2 Jy sample 115. In Figure ^ a tendency towards stronger clustering of the optically selected 
galaxies is seen on all scales. One has to bear in mind, that we compared two samples with only 
22% of the volume of the samples with 100/i~^Mpc depth. Hence it is not clear whether this is 
an effect only seen in our "local surrounding" . To do a reliable comparison of the morphology of 
optical and infrared selected galaxy distributions we have to wait for larger optical and infrared 
redshift surveys. 



7 Discussion and conclusions 

We analyzed morphological characteristics of the galaxy distribution given by the IRAS 1.2 Jy 
catalogue. The two subsamples (north and south) of this catalogue were studied individually with 
Minkowski functionals. 

Since the IRAS data have been obtained from a single instrument with uniform calibration, 
the two subsamples, which contain about the same number of galaxies, can be compared reliably. 

We have reported in detail our error estimates and have discussed tests on selection effects 
to assess the significance of our results. For reference purposes we used typical realizations of a 
stationary Poisson point process. 

Our results can be summarized as follows: The values of the Minkowski functionals for the 
southern part differ significantly from those for the northern part in the volume limited subsample 
with 100/i^^Mpc depth. This difference does not disappear on scales of 200/i~^Mpc as shown in 



Kerscher et al. (1997) 



Similar anisotropics in the angular distribution of IRAS galaxies around the northern and 
southern galactic poles have already been reported by Rowan-Robinson et al. (1986)| . However, 



the majority of previous IRAS catalogue studies have focused attention on the complete sample 
without addressing the distinction of its constituent parts. 

There is no reason to assume a distinguished position of the Milky Way galaxy; we therefore 
conclude that fluctuations in the global morphological characteristics of the IRAS sample ex- 
tend over length scales of at least 100/i~^Mpc. These fluctuations may originate from dynamical 
correlations in the matter dist ribution which arise during the cosmic e volution 



During the last few years, Coleman fc Pietronero (1992) , see also 3ylos Labini et al. (1997) 



in this volume, have advanced an interpretation of galaxy catalogue data in terms of a fractal 
support of the galaxy distribution. By its nature, a pure fractal indeed exhibits fluctuations in 
A'^-point distributions on all scales. Our results neither support nor contradict this interpretation, 
since the Minkowski functionals, as employed in the present paper, are global measures and are 
not designed to discriminate local structures of spatial patterns. 

Fluctuations occuring on scales up to lOOft-^^Mpc, at least, imply that cosmological simulations 
which enforce homogeneity on the scale of a few hundreds of /i~^Mpc and suppress fluctuations on 
larger scales by using periodic boundary conditions cannot reproduce the large-scale fluctuations 
indicated by the present analysis of the 1.2 Jy catalogue. This assertion is confirmed by our com- 



parison with IRAS mock catalogues drawn from simulations of 256/i ^Mpc box-length (Kerschei 



et al. 1997) 



It remains an open question, whether it is possible to explain these observed fluctuations in the 
clustering properties within large realizations of the standard model and with COBE normalized 
Gaussian initial density fields. 
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